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Abstract. Out of the four components of the electromagnetic field, Maxwell's theory only contains two physical degrees 
of freedom. However, in an expanding universe, consistently eliminating one of the "unphysical" states in the covariant 
(Gupta-Bleuler) formalism turns out to be difficult to realize. In this work we explore the possibility of quantization without 
subsidiary conditions. This implies that the theory would contain a third physical state. The presence of such a new (temporal) 
electromagnetic mode on cosmological scales is shown to generate an effective cosmological constant which can account for 
the accelerated expansion of the universe. This new polarization state is completely decoupled from charged matter, but can be 
excited gravitationally. In fact, primordial electromagnetic quantum fluctuations produced during electroweak scale inflation 
could naturally explain the presence of this mode and also the measured value of the cosmological constant. The theory is 
compatible with all the local gravity tests, it is free from classical or quantum instabilities and reduces to standard QED in 
the flat space-time limit. Thus we see that, not only the true nature of dark energy can be established without resorting to 
new physics, but also the value of the cosmological constant finds a natural explanation in the context of standard inflationary 
cosmology. Possible signals, discriminating this model from LCDM, are also discussed. 
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WHAT IS THE NATURE OF DARK ENERGY? 

Although the standard cosmological constant provides a simple and accurate phenomenological description of the 
present phase of accelerated expansion of the universe JD], it is well known that it suffers from an important naturalness 
problem. Indeed, the scale of the cosmological constant around 10~ 3 eV is many orders of magnitude smaller than the 
fundamental gravity scale given by the Planck constant which is around 10 19 GeV. 

Alternative models based on new physics have been proposed, however they are generally plagued by: 

• classical or quantum instabilities 

• fine tuning problems 

• inconsistencies with local gravity constraints 

Also, the possibility that the accelerated expansion could be signaling the breakdown of General Relativity on large 
scales has been also proposed, and several modified gravity theories have been considered 10, [H, although, in general, 
they also suffer from the same kind of problems described above. 

However, apart from gravity there is another long-range interaction in nature which could be relevant on cosmo- 
logical scales which is nothing but electromagnetism. It is usually assumed that due to the electric neutrality of the 
universe on large scales, the only relevant interaction in cosmology should be gravitation. However, the role played by 
electromagnetic astrophysical and cosmological fields is still far from clear, being the existence of cosmic magnetic 
fields the most evident example J3l. As a matter of fact, and as for any other interaction in nature, electromagnetism 
has only been tested in a limited range of energies or distances, which are roughly below 1.3 A.U. [5[] and, as a con- 
sequence, no experimental evidence is available for electromagnetic fields whose wavelenghts are comparable to the 
present Hubble radius. 

In this work, we consider (quantum) electromagnetic fields in an expanding universe. Unlike previous works, we 
will not follow the standard Coulomb gauge approach, but instead we will concentrate on the covariant quantization 
method. This method is found to exhibit certain difficulties when trying to impose the quantum Lorenz condition on 
cosmological scales [6]. In order to overcome these problems, and instead of introducing ghosts fields, we explore 
the possibility of quantizing without imposing any subsidiary condition [7, 8]. This implies that the electromagnetic 
field would contain an additional (scalar) polarization, similar to that generated in the massless limit of massive 
electrodynamics [9]. Interestingly, the energy density of quantum fluctuations of this new electromagnetic state 
generated during inflation gets frozen on cosmological scales, giving rise to an effective cosmological constant 0]. 



ELECTROMAGNETIC QUANTIZATION IN FLAT SPACE-TIME 



Let us start by briefly reviewing electromagnetic quantization in Minkowski space-time 1 10] since this will be useful 
in the rest of the work. The action of the theory reads: 

5=1/^-^^+4/) (1) 

where is a conserved current. This action is invariant under gauge tranformations — > + d^A with A an 
arbitrary function of space-time coordinates. At the classical level, this action gives rise to the well-known Maxwell's 
equations: 

d v F^ v =J^. (2) 

However, when trying to quantize the theory, several problems arise related to the redundancy in the description due 
to the gauge invariance. Thus, in particular, we find that it is not possible to construct a propagator for the A^ field 
and also that "unphysical" polarizations of the photon field are present. Two different approaches are usually followed 
in order to avoid these difficulties. In the first one, which is the basis of the Coulomb gauge quantization, the gauge 
invariance of the action (HJ is used to eliminate the "unphysical" degrees of freedom. With that purpose the (Lorenz) 
condition d^A^ = is imposed by means of a suitable gauge transformation. Thus, the equations of motion reduce to: 

□A M =J M . (3) 

The Lorenz condition does not fix completely the gauge freedom, still it is possible to perform residual gauge 
transformations A^ ^ A^ + 0, provided DO = 0. Using this residual symmetry and taking into account the form of 
equations (O, it is possible to eliminate one additional component of the A^ field in the asymptotically free regions 

(typically Aq) which means V -A — 0, so that finally the temporal and longitudinal photons are removed and we are left 
with the two transverse polarizations of the massless free photon, which are the only modes (with positive energies) 
which are quantized in this formalism. 

The second approach is the basis of the covariant (Gupta-Bleuler) and path-integral formalisms. The starting point 
is a modification of the action in (03, namely: 

S = Jd\ (™ V v + \{d^f +A^ . (4) 

This action is no longer invariant under general gauge transformations, but only under residual ones. The equations of 
motion obtained from this action now read: 



d v F^ v + ^{d v A v )=J^. (5) 

In order to recover Maxwell's equation, the Lorenz condition must be imposed so that the | term disappears. At 
the classical level this can be achieved by means of appropriate boundary conditions on the field. Indeed, taking the 
four-divergence of the above equation, we find: 

D(d v A v )=0 (6) 

where we have made use of current conservation. This means that the field d v A v evolves as a free scalar field, so that 
if it vanishes for large \t\, it will vanish at all times. At the quantum level, the Lorenz condition cannot be imposed 
as an operator identity, but only in the weak sense <9 V A V M \<j>) — 0, where (+) denotes the positive frequency part of 
the operator and \<j>) is a physical state. This condition is equivalent to imposing [ao(k) +aii(&)]|0) = 0, with ao and 
a the annihilation operators corresponding to temporal and longitudinal electromagnetic states. Thus, in the covariant 
formalism, the physical states contain the same number of temporal and longitudinal photons, so that their energy 
densities, having opposite signs, cancel each other. 

Thus we see that also in this case, the Lorenz condition seems to be essential in order to recover standard Maxwell's 
equations and get rid of the negative energy states. 



COVARIANT QUANTIZATION IN AN EXPANDING UNIVERSE 



So far we have only considered Maxwell's theory in fiat space-time, however when we move to a curved background, 
and in particular to an expanding universe, then consistently imposing the Lorenz condition in the covariant formalism 
turns out to be difficult to realize. Indeed, let us consider the curved space-time version of action ©: 



S= Id Xy/g 



-^ v ^ + i(V/) 2 +A/ 



(7) 



Now the modified Maxwell's equations read: 



(8) 



and taking again the four divergence, we get: 



□ (V V A V ) =0 



(9) 



We see that once again V V A V behaves as a scalar field which is decoupled from the conserved electromagnetic currents, 
but it is non-conformally coupled to gravity. This means that, unlike the flat space-time case, this field can be excited 
from quantum vacuum fluctuations by the expanding background in a completely analogous way to the inflaton 
fluctuations during inflation. Thus this poses the question of the validity of the Lorenz condition at all times. 

In order to illustrate this effect, we will present a toy example. Let us consider quantization in the absence of currents, 
in a spatially flat expanding background, whose metric is written in conformal time as ds 2 = a(r\) 2 (dr\ 2 — dx 2 ) with 
a{r\) — 2 + tanh(Tj/rjo) where rjo is constant. This metric possesses two asymptotically Minkowskian regions in the 
remote past and far future. We solve the coupled system of equations ([H) for the corresponding Fourier modes, which 
are defined as ^(r],x) = Jd 3 /b^(T7)e' K . Thus, for a given mode k, the stfn field is decomposed into temporal, 
longitudinal and transverse components. The corresponding equations read: 
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(10) 



with = a'/fl and k = \k\. We see that the transverse modes are decoupled from the background, whereas the 
temporal and longitudinal ones are non-trivially coupled to each other and to gravity. Let us prepare our system in 

an initial state \(j)} belonging to the physical Hilbert space, i.e. satisfying d v £/^ + ^\(j)) = in the initial flat region. 
Because of the expansion of the universe, the positive frequency modes in the in region with a given temporal or 
longitudinal polarization A will become a linear superposition of positive and negative frequency modes in the out 
region and with different polarizations A' (in this section we will work in the Feynman gauge | = — 1), thus we have: 



^X_(in) 
jik 



E 

A'=0,| 



n\ J^ 1 (out) n / j*\ (out) 



(11) 



or in terms of creation and annihilation operators: 



x=o,\\ 



(12) 



with A , A' = 0, || and where a%x' and j3;u' are the so-called Bogolyubov coefficients (see Hill for a detailed discussion), 
which are normalized in our case according to: 



E ( a Ap «A'p' ^PP' - Pxp h'p> *7pp') = TlXX' 

P,P'=0,\\ 



(13) 



n(k) 




FIGURE 1. Occupation numbers for temporal (continuous line) and longitudinal (dashed line) photons in the out region vs. k in 
rjQ 1 units. 



with r]ix' = diag{— 1,1) with X,X' — 0, ||. Notice that the normalization is different from the standard one fi"ill . 
because of the presence of negative norm states. 

Thus, the system will end up in a final state which no longer satisfies the weak Lorenz condition i.e. in the 

out region d v s^„Ji\^) ^ 0. This is shown in Fig. 1, where we have computed the final number of temporal and 
longitudinal photons n'£f(k) = Y.x lftu'(^)l 2 > starting from an initial vacuum state with = n'"(k) = 0. We see 

that, as commented above, in the final region n"" 1 (k) ^ n|"' (k) and the state no longer satisfies the Lorenz condition. 
Notice that the failure comes essentially from large scales, i.e. super-Hubble scales with kr\Q <C 1, since on small 
(sub-Hubble) scales with kr\Q 3> 1, the Lorenz condition is restored. This can be easily interpreted from the fact that 
on small scales the geometry can be considered as essentially Minkowskian. 

In order to overcome this problem, it is possible to impose a more stringent gauge-fixing condition. Indeed, we 
have shown above that in a space-time configuration with asymptotic fiat regions, an initial state satisfying the weak 
Lorenz condition does not necessarily satisfy it at a later time. However, it would be possible (see II 1210 to define the 
physical states |0) as those such that VnA^ + '\(j)) = 0, Vtj. Although this is a perfectly consistent solution, notice 
that the separation in positive and negative frequency parts depends on the space-time geometry and therefore, the 
determination of the physical states requires a previous knowledge of the future evolution of the universe. 

Another possible way out would be to modify the standard Gupta-Bleuler formalism by including ghosts fields as 
done in non-abelian gauge theories [13]. With that purpose, the action of the theory (0 can be modified by including 
the ghost term (see I14IP : 

S g = Jd 4 x^g^d^cd v c (14) 

where c are the complex scalar ghost fields. It is a well-known result lfl4l[l5ll that by choosing appropriate boundary 
conditions for the electromagnetic and ghosts Green's functions, it is possible to get (0|7]f v + T^ v \<j>) = 0, where T^ v 
and Tn V denote the contribution to the energy-momentum tensor from the % term in (Q and from the ghost term ( TT~4b 
respectively. Notice that a choice of boundary conditions in the Green's functions corresponds to a choice of vacuum 
state. Therefore, also in this case an a priori knowledge of the future behaviour of the universe expansion is required 
in order to determine the physical states. 

Although these are valid quanization procedures, in this work we follow a different approach in order to deal with the 
difficulties found in the Gupta-Bleuler formalism and we will explore the possibility of quantization in an expanding 
universe without imposing any subsidiary condition. 



QUANTIZATION WITHOUT THE LORENZ CONDITION 



In the previous section it has been shown that although the Lorenz gauge-fixing conditions can be formally imposed 
in the covariant formalism, this cannot be done in a straightforward way. These difficulties could be suggesting some 
more fundamental obstacle in the formulation of an electromagnetic gauge invariant theory in an expanding universe. 
As a matter of fact, electromagnetic models which break gauge invariance on cosmological scales have been widely 
considered in the context of generation of primordial magnetic fields (see, for instance, 1 16]). 

Let us then explore the possibility that the fundamental theory of electromagentism is not given by the gauge 
invariant action (Q~|i, but by the gauge non-invariant action (|7). Abandoning gauge invariance could, in principle, pose 
important problems for the viability of the theory, namely: 

• Modification of classical Maxwell's equations 

• Electric charge non-conservation 

• New unobserved photon polarizations 

• Negative norm (energy) states 

• Conflicts with QED phenomenology 

However, as we will show in the following, none of these problems is actually present for the theory in 

Since the fundamental electromagnetic theory is assumed to be non-invariant under arbitrary gauge transformations, 
then there is no need to impose the Lorenz constraint in the quantization procedure. Therefore, having removed one 
constraint, the theory contains one additional degree of freedom. Thus, the general solution for the modified equations 
^ can be written as: 

*p = ^ + *f 5 + + d »° ( 15 ) 

where srfjp with i = 1,2 are the two transverse modes of the massless photon, stfjf 1 is the new scalar state, which is the 
mode that would have been eliminated if we had imposed the Lorenz condition and, finally, dn is a purely residual 
gauge mode, which can be eliminated by means of a residual gauge transformation in the asymptotically free regions, 
in a completely analogous way to the elimination of the Aq component in the Coulomb quantization. The fact that 
Maxwell's electromagnetism could contain an additional scalar mode decoupled from electromagnetic currents, but 
with non-vanishing gravitational interactions, was already noticed in a different context in Jgl]. 

In order to quantize the free theory, we perform the mode expansion of the field with the corresponding creation and 
annihilation operators for the three physical states: 



J 1 1 i 



(16) 



A=l,2,s 

where the modes are required to be orthonormal with respect to the scalar product (see for instance H2]): 



S u ,S®(k-k'), A,A' = l,2,s (17) 



where dT.^ is the three-volume element of the Cauchy hypersurfaces. In a Robertson-Walker metric in conformal time, 
it reads c/E^ = a 4 (rj)(d 3 x, 0,0,0). The generalized conjugate momenta are defined as: 

Yl^ v = -(F» v -^g^VpAP) (18) 

Notice that the three modes can be chosen to have positive normalization. The equal-time commutation relations: 

[^(i7,f),^ v (77,x')] = [n°' i (Tj,J),n 0v (Tj ) j?/)] =o (19) 

and 

K(77,f),n 0v (rj,^)] =i^8®(x-x') (20) 



can be seen to imply the canonical commutation relations: 



a x (!),a£,(#)] =8 xv 5W(k-k>), A,A' = l,2,s (21) 

by means of the normalization condition in (T% . Notice that the sign of the commutators is positive for the three 
physical states, i.e. there are no negative norm states in the theory, which in turn implies that there are no negative 
energy states as we will see below in an explicit example. 

Since V^&fL 1 evolves as a minimally coupled scalar field, as shown in (0, on sub-Hubble scales (\kr]\ ^S> 1), we 

find that for arbitrary background evolution, \ oc a~ l , i.e. the field is suppressed by the universe expansion, 

thus effectively recovering the Lorenz condition on small scales. Notice that this is a consequence of the cosmological 
evolution, not being imposed as a boundary condition as in the flat space-time case. 

On the other hand, on super-Hubble scales (\kr\ | <C 1), \^n^jf \ — const, which, as shown in [7], implies that the 
field contributes as a cosmological constant in 0. Indeed, the energy-momentum tensor derived from (0 reads: 

T liv = -FiAaF v a + -guvFapF"! 3 



g^v 



(V a A a ) 2 +2A a V a (VpA?) 



4A (M V V) (V a A a ) 



(22) 



Notice that, for the scalar electromagnetic mode in the super-Hubble limit, the contributions involving F^ v vanish and 
only the piece proportional to S, is relevant. Thus, it can be easily seen that, since in this ca.se V^A 01 — constant , the 
energy-momentum tensor is just given by: 

^v = |^v(V a A a ) 2 (23) 

which is the energy-momentum tensor of a cosmological constant. Notice that, as seen in (0, the new scalar mode is a 
massless free field. This is one of the most relevant aspects of the present model in which, unlike existing dark energy 
theories based on scalar fields, dark energy can be generated without including any potential term or dimensional 
constant. 

Since, as shown above, the field amplitude remains frozen on super-Hubble scales and starts decaying once the mode 
enters the horizon in the radiation or matter eras, the effect of the % term in (0 is completely negligible on sub-Hubble 
scales, since the initial amplitude generated during inflation is very small as we will show below. Thus, below 1.3 AU, 
which is the largest distance scale at which electromagnetism has been tested [5], the modified Maxwell's equations 
(0 are physically indistinguishable from the flat space-time ones (0. 

Notice that in Minkowski space-time, the theory (0 is completely equivalent to standard QED. This is so because, 
although non-gauge invariant, the corresponding effective action is equivalent to the standard BRS invariant effective 
action of QED. Thus, the effective action for QED obtained from (0 by the standard gauge-fixing procedure reads: 

= / [dAp«^][^^ (24) 

where J£p is the Lagrangian density of charged fermions. The ^ term and the ghosts field appear in the Faddeev-Popov 
procedure when selecting an element of each gauge orbit. However, ghosts being decoupled from the electromagnetic 
currents can be integrated out in flat space-time, so that up to an irrelevant normalization constant we find: 



jw 



[dA^me 1 ^'^^^^ (25) 



which is nothing but the effective action coming from the gauge non-invariant theory (0 in flat space-time, in which 
no gauge-fixing procedure is required. 

To summarize, none of the above mentioned consistency problems for the theory in (0 arise, thus: 

• Ordinary Maxwell's equations are recovered on those small scales in which electromagnetism has been tested. 

• Electric charge is conserved since only the gauge electromagnetic sector is modified but not the sector of charged 
particles which preserves its gauge symmetry. 

• The new state only couples gravitationally and evades laboratory detection. 



• The new state has positive norm (energy). 

• The effective action is completely equivalent to standard QED in the fiat space-time limit. This guarantees, not 
only that the standard phenomenology is recovered, but also that no new interaction terms will appear in the 
renormalization procedure. 



COSMOLOGICAL EVOLUTION 

Let us now consider the cosmological evolution of this new electromagnetic mode. For that purpose, we will con- 
sider an homogeneous electromagnetic field with = (Ao(t),A(t)). The corresponding equations motion in a fiat 
Robertson- Walker background in cosmological time t read: 

A + 3HA + 3HA = 

A + HA = (26) 

Notice that conformal time components are related to components in cosmological time as stf^ = (oAq,A). The effects 
of the high electric conductivity of the universe can be taken into account by including the corresponding current term 
on the right hand side of Maxwell's equations, however due to the strict electric neutrality of the universe on large 
scales = 0. 

In the case in which the scale factor behaves as a simple power law with H = p/t, the solutions for the above 
equations are: 

A (t) = A+t+A r 3 P 

A(t) = A + t l ~ p +A~ (27) 

Thus we see that the temporal component grows faster than the spatial one and therefore, at late times, on cosmological 
scales we can ignore the spatial contribution and the new scalar state is essentially given by the electric potential. Notice 
also that from the temporal equation we get: 

j t (V^) = j t (A + 3HA Q )=0 (28) 

irrespective of the background evolution, i.e. on large scales the scalar mode V^A^ 1 is strictly constant as commented 
before. We can also compute the contributions from the temporal and spatial components to the energy density from 
d22l i, thus we get: 

A • 2 
Pa = ^{Ao + IHAq) = const. 

PA = i(^) 2 ~^ 4 (29) 

We see that also, as commented before, the temporal component behaves as a cosmological constant, whereas the 
contibution from the spatial part decays as radiation and therefore does not affect the universe isotropy on large scales. 

In Fig. 2 we show the cosmological evolution of the electric potential Aq. We see that the field is constant during 
inflation, it grows linearly in time in the matter and radiation eras and becomes also constant when the eletromagnetic 
dark energy starts dominating. Notice that the present value of Aq will be determined by the initial value of the field 
generated during inflation from quantum fluctuations. This amplitude can be explicitly calculated as shown in next 
section. 



QUANTUM FLUCTUATIONS DURING INFLATION 



Let us consider the quantization in an inflationary de Sitter space-time with fl(fj) = — l/(H]T]). Here we will take 
£, = 1/3, although similar results can be obtained for any | > 0. The explicit solution for the normalized scalar state 
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where Ei(jc) = /"e is the exponential integral function. Using this solution, we find: 

v - _ a ~ 2 ^) 3 -&!,+& 

* (2^)3/2 ^21 2 fc^ 

so that the field is suppressed in the sub-Hubble limit as V^jz^y ~ £?((krj)~ 2 ). 

On the other hand, from the energy density given by p A = T° Q , we obtain in the sub-Hubble limit the corresponding 
Hamiltonian, which is given by: 



H 



(32) 



We see that the theory does not contain negative energy states (ghosts). 

Finally, from ( f3Qb it is possible to obtain the dispersion of the effective cosmological constant during inflation: 

dk , 



<0|(V^) 2 |0)= /y^(*) 



with P A (k) = 47T/t 3 |V M ^ M ' 1 | 2 . In the super-Hubble limit, we obtain for the power-spectrum: 



P A {k) 
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(33) 



(34) 



in agreement with @|. Notice that this result implies that p A ~ {Hi) 4 . The measured value of the cosmological constant 
then requires Hj ~ 10~ 3 eV, which corresponds to an inflationary scale Mj ~ 1 TeV. Thus we see that the cosmological 
constant scale can be naturally explained in terms of physics at the electroweak scale. 



PERTURBATIONS AND CONSISTENCY 



Despite the fact that the background evolution in the present case is the same as in ACDM, the evolution of metric 
perturbations could be different, thus offering an observational way of discriminating between the two models. With 
this purpose, we have calculated the evolution of metric, matter density and electromagnetic perturbations [17]. The 
propagation speeds of scalar, vector and tensor perturbations are found to be real and equal to the speed of light, 
so that the theory is classically stable. We have also shown before that the theory does not contain ghosts and it is 
therefore stable at the quantum level. On the other hand, using the explicit expressions in 111 811 for the vector-tensor 
theory of gravity corresponding to the action in (0, it is possible to see that all the parametrized post-Newtonian (PPN) 
parameters agree with those of General Relativity, i.e. the theory is compatible with all the local gravity constraints for 
any value of the homogeneous background electromagnetic field IT7L fl9fl . 

Concerning the evolution of scalar perturbations, we find that the only relevant deviations with respect to ACDM 
appear on large scales k ~ Hq and that they depend on the primordial spectrum of electromagnetic fluctuations. 
However, the effects on the CMB temperature and matter power spectra are compatible with observations except 
for very large primordial fluctuations. In addition, the different evolution of the scalar potential <f><. with respect to the 
ACDM model gives rise to a possible discriminating contribution to the late-time integrated Sachs-Wolfe effect |20ll . 

CONCLUSIONS 

In this work we have presented a consistent quantization procedure for the electromagnetic interaction with three 
polarization states. The energy density of the new scalar mode on cosmological scales is shown to behave as an 
effective cosmological constant, whose value is determined by the amplitude of quantum fluctuations generated during 
inflation. As a matter of fact, the measured value of the cosmological constant is naturally explained provided inflation 
took place at the electroweak scale. The model is free from classical or quantum instablities and is consistent with 
all the local gravity constraints. On the other hand, it is also compatible with observations from CMB and large scale 
structure and contains the same number of free parameters as ACDM. Unlike dark energy models based on scalar fields, 
acceleration in this model arises from the kinetic term of the new electromagnetic mode, without the introduction of 
unnatural dimensional parameters or potential terms. The results presented in this work show that, unlike previous 
proposals, the nature of dark energy can be established without resorting to new physics. 
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